Beam quality of a non-ideal atom laser 
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We study the propagation of a non-interacting atom laser distorted by the strong lensing effect of 
the Bose-Einstein Condensate (BEC) from which it is outcoupled. We observe a transverse structure 
containing caustics that vary with the density within the residing BEC. Using WKB approximation, 
Fresnel-Kirchhoff integral formalism and ABCD matrices, we are able to describe analytically the 
atom laser propagation. This allows us to characterize the quality of the non-ideal atom laser beam 
by a generalized M 2 factor defined in analogy to photon lasers. Finally we measure this quality 
factor for different lensing effects. 
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Optical lasers have had an enormous impact on 
science and technology, due to their high brightness 
and coherence. The high spatial quality of the beam 
and the little spread when propagating in the far-field 
enable applications ranging from the focusing onto 
tiny spots and optical lithography [l| to collimation 
over astronomic distances [2|. In atomic physics, Bose- 
Einstein condensates (BEC) of trapped atoms Q are an 
atomic equivalent to photons stored in a single mode of 
an optical cavity, from which a coherent matter wave 
(atom laser) can be extracted EL El- The possibility of 
creating continuous atom laser |f]J promises spectacular 
improvements in future applications 0, 
where perfect collimation or strong focusing [121 \1' 
are of prior importance. Nevertheless these properties 
depend drastically on whether the diffraction limit can 
be achieved. Thus, characterizing the deviation from this 
limit is, as for optical lasers 15], of crucial importance. 
For example, thermal lensing effects in optical laser 
cavities, which cause significant decollimation, can also 
induce aberrations that degrade the transverse profile. 
In atom optics, a trapped BEC weakly interacting with 
the outcoupled atom-laser beam acts as an effective 
thin-lens which leads to the divergence of the atom laser 
[lfj| without affecting the diffraction limit. When the 
lensing effect increases, dramatic degradations of the 
beam are predicted 17], with the apparition of caustics 
on the edge of the beam. 

In order to quantitatively qualify the atom-laser beam 
quality, it is tempting to take advantage of the methods 
developed in optics to deal with non- ideal laser beams i.e. 
above the diffraction limit. Following the initial work of 
Siegman [18( who introduced the quality factor M 2 which 





FIG. 1: Absorption images of a non-ideal atom laser, cor- 
responding to density integration along the elongated axis 
of the BEC. The figures correspond to different height of 
RF-outcoupler detunings with respect to the bottom of the 
BEC: (a) 0.37 fim (b) 2.22 ^m (c) 3.55 fim. The graph above 
shows the RF-outcoupler (dashed line) and the BEC slice 
(red) which is crossed by the atom laser. This results in the 
observation of caustics. The field of view is 350 /im x 1200 /im 
for each image. 



is proportional to the space-beam-width (divergence x 
size) product at the waist, it is natural to extend its 
definition to atom optics as 



M 2 

Ax Ak x = — , 



(1) 
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where Ax and Ak x = Ap x /h characterize respectively 
the size and the divergence along x (Ap x is the width 
of the momentum distribution). Equation plays the 
same role as the Heisenberg dispersion relation: it ex- 
presses how many times the beam deviates from the 
diffraction limit. 

In this letter, we experimentally and theoretically 
study the quality factor M 2 of a non-ideal, non- 
interacting atom-laser beam. First, we present our exper- 
imental investigation of the structures that appear in the 
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transverse profile 0] . We show that they are induced by 
the strong lensing effect due to the interactions between 
the trapped BEC and the outcoupled beam. Then, us- 
ing an approach based on the WKB approximation and 
the Fresnel-Kirchhoff integral formalism, we are able to 
calculate analytical profiles which agree with our exper- 
imental observations. This allows us to generalize con- 
cepts introduced in for photon laser and to calculate 
the quality factor M 2 . This parameter can then be used 
in combination with the paraxial ABCD matrices to 
describe the propagation of the non-ideal beam via the 
evolution of the rms width. Finally, we present a study 
of the M 2 quality factor as a function of the thickness of 
the BEC-induced output lens. 

Our experiment produces atom lasers obtained by ra- 
dio frequency (RF) outcoupling from a BEC |5j, |2lJ . The 
experimental setup for creating condensates of 87 Rb is de- 
scribed in detail in • Briefly, a Zeeman-slowed atomic 
beam loads a magneto-optical trap in a glass cell. About 
2 x 10 s atoms are transferred in the \F,mp) = |1, — 1) 
state to a Ioffc-Pritchard magnetic trap , which is sub- 
sequently compressed to oscillation frequencies of u) y — 
2ir x 8 Hz and uj x>z = 2tt x 330 Hz in the dipole and 
quadrupole directions respectively. A 25 s RF-induced 
evaporative cooling ramp results in a pure condensate of 
N = f0 6 atoms, cigar-shaped along the y axis. 

The atom laser is extracted from the BEC by applying 
a RF field a few kHz above the bottom of the trap, in 
order to couple the trapped state to the weakly anti- 
trapped state |1,0). The extracted atom laser beam falls 
under the effect of both gravity —mgz and second order 
Zeeman effect V = —muj 2 (x 2 + z 2 )/2 with u = 2tt x 
20 Hz (see Fig. EH)- The RF-outcoupler amplitude is 
weak enough to avoid perturbation of the condensate so 
that the laser dynamics is quasi-stationary [2l| and the 
resulting atom flux is low enough to avoid interactions 
within the propagating beam. Since the BEC is displaced 
vertically by the gravitational sag, the value of the RF- 
outcoupler frequency i/rp defines the height where the 
laser is extracted 5]- After f ms of operation, the fields 
are switched off and absorption imaging is taken after 1 
ms of free fall with a measured spatial resolution of 6 
fim. The line of sight is along the weak y axis so that 
we observe the transverse profile of the atom laser in the 
(z, x) plane. 

Typical images are shown in figure^ Transverse struc- 
tures, similar to the predictions in |l7j . are clearly visible 
in figures QJd and^. The laser beam quality degrades as 
the RF-outcoupler is higher in the BEC (i.e. the laser 
beam crosses more condensate) , supporting the interpre- 
tation that this effect is due to the strong repulsive inter- 
action between the BEC and the laser. This effect can be 
understood with a semi-classical picture. The mean-field 
interaction results in an inverted harmonic potential of 
frequencies u>i (in the directions i = x,y, z) which, in the 
Thomas-Fermi regime, are fixed by the magnetic confinc- 
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FIG. 2: (a) Principle of the RF-outcoupler : the radio- 
frequency vrf = (-Ebec — Eiasei)/h selects the initial position 
of the extracted atom-laser beam. The laser is then subjected 
to the condensate mean- field potential (Tint, to the quadratic 
Zeeman effect V and to gravity —mgz. For the sake of clar- 
ity, V has been exaggerated on the graph, (b) Representation 
of the two-dimensional theoretical treatment: (I) inside the 
condensate, phase integral along atomic paths determines the 
laser wavefront at the BEC output (WKB approximation). 
(II) A Fresnel-Kirchhoff integral on the contour F is used to 
calculate the stationary laser wavefunction at any point below 
the condensate. (Ill) As soon as the beam enters the parax- 
ial regime, we calculate the M 2 quality factor and use ABCD 
matrix formalism. 



ment (lill^ . The interaction potential expels the atoms 
transversally, as illustrated in Fig. [2t> - Because of the fi- 
nite size of the condensate, the trajectories initially at 
the center of the beam experience more mean-field repul- 
sion than the ones initially at the border. This results 
in accumulation of trajectories at the edge of the atom 
laser beam |l7| . in a similar manner to caustics in op- 
tics. This picture enables a clear physical understanding 
of the behaviour observed in figure ^ if i^p is chosen 
so that extraction is located at the bottom of the BEC 
(Fig. [IJi), the lensing effect is negligible and one gets a 
collimated beam. As the RF outcoupler moves upwards 
(Fig. ^ and ^0, a thicker part of the condensate acts 
on the laser and defocusing, then caustics appear. We 
verified that when sufficiently decreasing the transverse 
confinement of the trapped BEC, i.e. making the inter- 
action with the outcoupled atoms negligible, the atom 
laser is collimated at any RF value |2(j. 

In order to describe quantitatively the details of the 
profiles of the non-ideal atom laser one could solve nu- 
merically the Gross-Pitaevskii equation (CPE) [tR H^|. 
As we show here, another approach is possible, using ap- 
proximations initially developed in the context of photon 
optics and extended to atom optics. These approxima- 
tions allow calculation of the atom-laser propagation to- 
gether with the characterization of its rms width evolu- 
tion by means of the quality factor M 2 , in combination 
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with ABCD matrices. Note however that these matrices 
can only be used in the paraxial regime, i.e. when the 
transverse kinetic energy is smaller than the longitudinal 
one ^(|. This condition does not hold in the vicinity of 
the BEC, and thus we split the atom laser evolution into 
three steps (Fig. |2t>) where we use different formalisms 
in close analogy with optics: (I) WKB inside the con- 
densate (eikonal), (II) Fresnel-Kirchhoff integral outside 
the condensate and, (III) paraxial ABCD matrices after 
sufficient height of fall. In all the following, we do not 
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FIG. 3: Beam profile after 600 of propagation for a out- 
coupler height of 3.55 fira: (a) calculated in the central plane 
y = 0; (b) resulting from the integration along the line of 
sight y; (c) obtained experimentally. Deviations between the 
curves (b) and (c) can be attributed to imperfections in the 
imaging process and bias fluctuations. Note that for all three 
profiles, the overall shape is preserved and the calculated rms 
width in the central plane is in agreement with the measured 
one within experimental errors. 

include interactions within the atom laser since they are 
negligible for the very dilute beam considered in this let- 
ter. In addition, the condensate is elongated along the y 
axis, so that the forces along this direction are negligible 
and we consider only the dynamics in the (z, x) plane. 

First, we consider that the atoms extracted from the 
condensate by the RF-outcoupler start at zero velocity 
from ro. In region (I) the total potential (resulting from 
gravity and interactions with the trapped BEC) is cylin- 
drically symmetric along the y axis. The trajectories 
are thus straight lines (see Fig. The beam profile 

?/>(ri) at the BEC border is obtained thanks to the WKB 
approximation, by integrating the phase along classical 
paths of duration To 

^( ri ) a / • i !n T etJr ° k(r) ' dr ^BEc(r ), (2) 
y/smh (2lj z to) 

where "0bec is the condensate wavefunction and the pref- 
actor ensures the conservation of the flux. 

The atom laser being in a quasi-stationary state, the 
wavefunction satisfies in region (II) the time- independent 
Schrodinger equation at a given energy E in the potential 
V— mgz. Since this equation is analogous to the Helmoltz 
equation j^, in optics, the Fresnel-Kirchhoff integral 
formalism can be generalized to atom optics as |29j 

V>(r) oc j dh • [G s ViV(ri) - V(ri)ViGjs] ■ (3) 

We take the contour T along the condensate border and 
close it at infinity (see Fig. ^jp ). The time- independent 



Green's function Gs(r, ri) is analytically evaluated from 
the time-domain Fourier transform of the Feynman prop- 
agator K(r, ri,r) [23,|33 calculated by means of the van 
Vleck formula [31|. Using ip( r i) from equation (J2J, the 
wavefunction ^(r) is then known at any location r. We 
verified the excellent agreement of this model with a nu- 
merical integration of the GPE including intra-laser in- 
teractions, thus confirming that they remain negligible 
throughout propagation. 

The method using a Kirchhoff integral is demanded 
only for the early stages of the propagation. As soon 
as the paraxial approximation becomes valid, the ABCD 
matrix formalism can be used to describe the propaga- 
tion of the beam. An example of a profile calculated in 
the central plane (y = 0) is presented in figure^. In 
figure |3t> we add the profiles of all y planes taking into 
account the measured resolution of the imaging system. 
In figure we present the corresponding experimental 
profile. The overall shape is in good agreement with the- 
ory, and the differences can be explained by imperfec- 
tions in the imaging process and bias fluctuations during 
outcoupling. The rms size of the three profiles, which 
depend very smoothly on the details of the structure, are 
the same within experimental uncertainties. Thus, here- 
after, as in the case of propagation of optical laser using 
the M 2 parameter, we only consider the evolution of the 
rms size of the atom laser. 

To calculate the beam width change in the paraxial 
regime, we define, following 18j, a generalized complex 
radius of curvature 

f iM 2 

where a x is the rms width of the density profile, a 
reduced variable which describes the time evolution of 
the beam such that £ = corresponds to the position of 
the waist. Equation involves an invariant coefficient, 
the beam-quality factor M 2 , as defined in Eq. 0J. This 
coefficient, as well as the effective curvature C(£), can be 
extracted from the wavefront in the paraxial domain, as 
explained in |32| . In optics, the generalized complex ra- 
dius obeys the same ABCD propagation rules as does a 
Gaussian beam of the same real beam size, if the wave- 
length A is changed to M 2 A |2jj. Similarly, the complex 
radius follows here the ABCD law for matter- waves 
|l6|. and we obtain the rms width 

2t** 2 ,2^ fM 2 h\ 2 sinh 2 (0 

where a x o is taken at the waist. 

This generalized Rayleigh formula allows us to measure 
M 2 . In the inset of figure the evolution of the trans- 
verse rms width a x versus propagation £, taken from ex- 
perimental images, is compared to the one given by Eq. 
|J5J) , where a x o is calculated with our model. For a chosen 
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RF-outcoupler position, we fit the variation of the width 
with a single free parameter M 2 . We then plot the mea- 
sured value of M 2 vs the ouput coupler height (Fig. 0}, 
and we find good agreement with theory. 
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FIG. 4: M 2 quality factor vs RF-outcoupler distance from the 
bottom of the BEC: theory (solid line), experimental points 
(circles). The two diamonds represent the M 2 for the two 
non-ideal atom lasers shown in figure 03 and 0;. The RF- 
outcoupler position is calibrated by the number of outcoupled 
atoms. Inset: typical fit of the laser rms size with the gener- 
alized Rayleigh formula (Eq. |SJ for RF-outcoupler position 
3.55 fim. 

In conclusion, we have characterized the transverse 
profile of an atom laser. We demonstrated that, in our 
case, lensing effect when crossing the condensate is a crit- 
ical contributor to the observed degradation of the beam. 
We showed that the beam-quality factor M 2 , initially in- 
troduced by Siegman for photon laser, is a fruitful 
concept for describing the propagation of an atom laser 
beam with ABCD matrices, as well as for characteriz- 
ing how far an atom laser deviates from the diffraction 
limit. For instance, it determines the minimal focusing 
size that can be achieved with atomic lenses provided 
that interactions in the laser remain negligible [lil . 
This is of essential importance in view of future appli- 
cations of coherent matter-waves as, for example, when 
coupling atom lasers onto guiding structures of atomic 
chips |34| . In addition, if interactions within atom laser 
become non negligible, a further treatement could be de- 
veloped in analogy with the work of 35] for non-linear 
optics. 
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